ERROR DISTRIBUTIONS FOR RANDOM GRID 
APPROXIMATIONS OF MULTIDIMENSIONAL 
STOCHASTIC INTEGRALS 

CARL LINDBERG 1 , HOLGER ROOTZEN 2 

Abstract. This paper proves joint convergence of the approximation 
error for several stochastic integrals with respect to local Brownian 
semimartingales, for non-equidistant and random grids. The conditions 
needed for convergence are that the Lebesgue integrals of the integrands 
tend uniformly to zero and that the squared variation and covariation 
processes converge. The paper also provides tools which simplify check- 
ing these conditions and which extend the range for the results. These 
results are used to prove an explicit limit theorem for random grid ap- 
proximations of integrals based on solutions of multidimensional SDE-s, 
and to find ways to "design" and optimize the distribution of the ap- 
proximation error. As examples we briefly discuss strategies for discrete 
option hedging. 



1. Introduction 

The error in numerical approximations of stochastic integrals is a random 
variable, or, if one also is interested in the "time" development of the error, 
a stochastic process. Hence the most precise evaluation of the error which is 
possible to obtain is to derive the distribution of the error. The prototype 
example is the Euler method for the stochastic integral J" * f(B(s), s)dB(s), 
for a Brownian motion B. The Euler method is to approximate the integrand 
with a stepfunction which is constant between the "evaluation times" (or in 
finance terminology "intervention times") of the grid i/n;i = 0, 1, . . . . This 
leads to the appro ximation L f or ]„dB(s), with r/ n (i) = i/n on the intervals 



i/n, (i + l)/n). In Rootzenl (1980) it is shown that the approximation error 



Date: August 26, 2011. 

2010 Mathematics Subject Classification. Primary 60F05, 60H05, 91G20; Secondary 
60G44, 60H35. 

Key words and phrases. Continuous local martingale, Brownian motion, random grid, 
error distribution, Euler method, weak convergence, optimality. 

1 carl. lindberg@alumni.chalmers.se Work supported by the Swedish Foundation for 
Strategic Research through the Gothenburg Mathematical Modelling Centre. 

2 hrootzen@chalmers.se Work supported by the Swedish Foundation for Strategic Re- 
search through the Gothenburg Mathematical Modelling Centre. 

1 



2 



CARL LINDBERG, HOLGER ROOTZEN 



9f(x,y) 



jjn _ n i/2 jt^j? _ j Q iq^dB^ converges stably in distribution, 

CT^ s i=y* f'(B(s),s)dW(s), 

where W is a Brownian motion independent of B and f'(x, y) 
and where Renyi's quite useful concept of stable convergence means that 
U n converges jointly with any sequence which converges in probability. 

The intuition behind this result is that "the small wiggles of a Brownian 
path are asymptotically independent of the global behaviour of the path" . 
The result has seen much further development, in particular to the error 
in numerical solution schemes for SDE's, and has recently found significant 
application in measuring the risks associated with discrete hedging. A brief 
overview of some of this literature is given below. 

The present paper generalizes this result in three ways: to joint con- 
vergence of the approximation error for several stochastic integrals, to lo- 
cal Brownian semimartingales instead of Brownian motions, and to non- 
equidistant and random evaluation times. The tools which help us quantify 
the intuition given above is Girsanov's theorem which shows how a mul- 
tidimensional Brownian motion is affected by a change of measure, and 
Levy's characterization of a multidimensional Brownian motion in terms of 
its square variation processes. 

The conditions needed for convergence apply more generally than to ap- 
proximation schemes. They are that the Lebesgue integrals of the integrands 
tend uniformly to zero in probability and that the square variation and co- 
variation processes converge in probability. We additionally provide tools 
which simplify checking these conditions and which extend the range of the 
results. Further we apply these results to prove an explicit limit theorem for 
approximations of integrals based on solutions of multidimensional SDE-s. 

One center of interest for this paper is the possibility to improve ap- 
proximation by using variable and random grids. In particular we study 
approximation schemes where the evaluation times i/n are replaced by time 
points given by the recursion Tq = and 

n n j ^ 

r k+1 -r k + ne(j ny 

for a positive adapted process 9{t). We also study how the function 9 can be 
chosen to design the approximation error so that it has desirable properties. 
For example, these could be homogeneous evolution of risk, or to make the 
approximation error have minimal standard deviation. 

A main motivation for writing this paper was to provide tools to study 
discrete hedging which uses random intervention times. We exemplify these 
possibilities by using the general results to exhibit a " no bad days" strategy 
and a minimum standard deviation strategy for the Black-Scholes model. 

Weak convergence theory for approximations of stochas tic integrals an d 
solutions to stochastic differential equations is developed in lRootzin 
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Kurtz and Protter ljl991«J ). lKurtz and Protterl (jl991bh . and lKurtz and Protterl 
(|l99ld ) and, in particular, an exten s ive st udy of the Euler method for SDE- 
s is provided by Jacod and Protter ( 19981 ). This theory has been used and 
extended to solve and analyze various aspects of approximation and hedg- 
ing error problems in M at hematical Finan ce. As examples we mention 

Duffie and Protterl (ll992T). iBertsimas et all (l200Ch. lHavashi and Mvkland 



( 200 51) .iTankov and Voltchkoval (|200fll ) . iBroden and Wiktorssonl (j20ld ). and 
FukasawaT(2011). A Malliav in calculus approach to discrete hedging is used 
in I Gob et and Temam ( 200ll ) and in a number of papers, w hich also consider 
variab le but deterministic grids, by Geiss and coworkers, see lGeiss and Toi vola 



(2009) and the references therein. The main theoretical tool of lHavashi and Mvklandl 



(|2005l ) is related t o our gene r al res ult, as discussed further below. The quite 
interesting paper iFukasawal l|201ll ) also studies random grid grid approxi- 
mations, for one-dimensional processes. The setting of Fukasawa's paper is 
more or less in the middle between our theorems 12. 21 and 13.31 The conditions 
used by Fukasawa are rather different from ours, and there doesn't seem to 
be any simple relations between the his results and ours. 

Now a brief overview of the paper. The next section, Section [21 contains 
the basic general theorem on multidimensional convergence for stochastic 
integrals with respect to local multidimensional Brownian semimartingales, 
and the tools to check conditions and extend the result. In Section [3] we give 
the explicit result for random grid approximations of stochastic integrals 
based on the solution of a multidimensional SDE. Section 0] investigates 
ways to design and optimize approximation errors, and in Section [5] this is 
applied to discrete financial hedging. 



2. General results 

This section contains two main results. The first one gives a means to 
establish multidimensional convergence of the distribution of stochastic in- 
tegrals with more and more rapidly varying integrands, and the second one 
shows how convergence of integrals with simple integrands can be extended 
to more general integrands. In addition, Lemma 12.81 provides tools to check 
the assumptions of the theorems. Our main aim is the error in approxi- 
mations of stochastic integrals, but the reults may in fact also have more 
general use. 

Let Q = C(M+,K d ) be the space of continuous Revalued functions de- 
fined on R + , define B t = {-Bt}i<«<d by Bt{<jj) = u>(t), let P be the probability 
measure which makes B a Brownian motion starting at 0, and let Tt be the 
completion of the cr-algebra generated by {B s ;0 < s < t}. Further write 
T for the smallest sigma algebra which contains all the Tt- Until further 
notice is given all random variables we consider are defined on the filtered 
probability space (0, {Tt)-, T, P). Weak convergence will be for random vari- 
ables (or "processes") with values in C([0,T],~R K ), the space of continuous 
X-dimensional functions defined on the the time interval [0, T], and with 
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respect to the uniform metric. Usually the dimension K of the processes 
will be clear from the context, and then we for brevity write C[0, T], instead 
of C([0,T],W K ), and just write for weak convergence. 

Weak convergence is stable (or " Renyi-stable" ) if it holds on any subset 
of J 7 , and the convergence is mixing ( or " Renyi- mixing" ) if in addition the 
limit is the same on any subset. In the present setting this is specified by 
the following definition. 

Definition 2.1. Let (X n ) n >i be a sequence of random variables defined 
on the same probability space (f^J 7 , P) and with values in C[0,T]. Then 
(X n ) n >i converges stably to X, written X n => s X, if E[Uf(X n )] con- 
verges for any bounded continuous function f : C[0, T] — > M. and any bounded 
measurable random variable U defined on (£l,J~, P). If in addition 

(1) ]imE[Uf(X n )]=E[U]E[f(X)}, 

n 

then (X n ) n >\ converges mixing to X, written X n => TO X. 

It is straightforward to see that to establish stable or mixing convergence 
it is enough to prove convergence of E[Uf(X n )] for strictly positive U with 
EU = 1. 

Let X* be a C[0, T]-valued random variable, defined on an extension 
(J7, J 7 , P) of the original probability space , such that X* has the same distri- 
bution as X and is independent of any J-"-measurable random variable. Then 
(see e.g. lAldous and Eaglesonl (Il978l ^ X n =^ m X iff{Y n ,X n ) => (Y,X*) 



for any sequence of random variables Y n — k, Y which converges in proba- 
bility iff X n => X with respect to ¥(-\A) for any set A with F(A) > 0. (In 
the middle statement, convergence is with respect to the product topology.) 
The same equivalencies hold for stable convergence, if the requirement "X* 
is independent of any J 7 - measurable random variable" is removed. 

Let X = (Xj,j = l,...d) be a continuous ci-dimensional local semi- 
martingale defined on the space (Q, (J^), P) by 



(2) 



d i-t r-t 

x A t ) = Y\ / G j)k {s)dB k (s) + / aj (s)ds, 
k=1 Jo Jo 



with Gj t k and aj adapted, and with J* Q T k ds < oo and ajds < oo a.s. 
for all j, k. Further let {Hfj} = {Hff 1 < i, j < d} he & d x d -dimensional 

array of J^-adapted processes such that f^(H^j) 2 dt < oo a.s. for each i,j, 
and write 

(3) {H?j.Xj} = {H^.X f ,l<i,j<d} 

= { [ H^(s)dX 3 (s); 1 < i,j < d} < t < T . 
Jo 

Thus {Hfa ■ Xj} takes values in C([0, T], R dxd ). Let -t p denote convergence 
in probability. 
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Theorem 2.2. Suppose that {Hf^} satisfies 
(4) sup | / H™jds\ -> p 0, n — > oo, 



o<t<T Jo 

and that 



n — > oo, 



(5)/ H™jGj j kHi' rn G m ^kds -^■ p f (:j HijGj^Hi rn G m ^Piji m ds, 
Jo 

for i,j,k,l,m = l,...,d, and for some correlation matrix processes p k = 
(Pij j m );^ = l,"^ an d positive process {Hij : 1 < i,j < d}. Let a k be an 
arbitrary roo t of p k such that a k (a k )' = p k for all t G [0,T]. Then, for X 
given by ([2]) ; 

d d 
k=l l,m=l 

as n — > oo, where W = (W^ m ^', 1 < l,m,k < d) is a d x d x d-dimensional 
Brownian motion which is independent of J-. 

This result simplifies in the special case when X is just a Brownian motion 
B, see the following corolla ry. The corollary is close to Theorem A.l of 
Havashi and Mvklandl (120051 ). Differences are that the corollary makes the 



basic condition ^ explicit, gives a more detailed description of the limit 
distribution, and has the more powerful conclusion of stable convergence. 

Corollary 2.3. Suppose that (j3|) is satisfied and that 

(6) f HfjHZ jds -> p [ H id H k jp> ik ds, n oo, 
Jo ' Jo 

as n —7- oo ; for some correlation matrix processes p 3 = a 3 (a 3 )', where 
i,j,k = 1, . . . d and positive processes {Hi j : 1 < i,j < d}, and < t < T. 
Then 

d 

(7) {u-r'h} 

k=l 

as n — > oo, where W = {Wjk\ 1 ^ k, j < d) is a dx d-dimensional Brownian 
motion which is independent of T . 

The following lemma plays in important role in the proofs. 

Lemma 2.4. Suppose that rj(t) and H n (t) are real-valued random processes 
with f r/(t) 2 dt < oo a.s. and with lim sup n _ >00 f Q H n (t) 2 dt < oo a.s. for 
some positive constant S < oo. Suppose further that 



sup | / H n ds\ — > p 0, // — x . 
o<t<s Jo 



6 



CARL LINDBERG, HOLGER ROOTZEN 



Then 

(8) sup I / H n r]ds\ —?- p 0, n — > oo. 

o<t<s Jo 

Proof of Lemma \2Ji\ Suppose first that there exists a sequence {%} of pro- 
cesses such that 



o 



s 

(rj(t) — 7]k(t)) 2 dt — > p as k 



oo, 



sup I / H n n k {s)ds\ — > p as n — > oo, for each k. 
Kt<s Jo 



o<t<s Jo 

Then, by the Cauchy-Schwarz inequality 



limsup sup | / H n t]ds\ < limsup sup | / H n ri k ds\ 

n 0<t<S Jo n 0<i<5 Jo 

+ limsup sup | / H n (r] — r]k)ds\ 
n 0<t<S Jo 



< O + Wlimsupy (H n ) 2 dtJ j (rj - f] k ) 2 dt 



which tends to as k — > oo, so that ([8]) holds. 

Thus the lemma follows if there exist a sequence {ijk} which satisfies 
the two requirements above. Now, e.g. by convolving rj with a sufficiently 
narrow normal density it can be seen that there exists a continuous process 
fjk, jointly measurable in t and u, such that F(f^ (n(t) — fjk(t)) 2 dt > 1/k) < 
1/k. Furthermore, with 1a denoting the indicator function of a set A, for 
Vk,m(t) = E!=o' Vk(iS/m)l te[iS/m ^ i+1)s/m) it follows that 
s 

(.Vk(t)-fjk,m(t)) 2 dt^ a , s .Q as m->oo, 

o 

and thus, choosing m k suitably, rj k = Y/^o Vk{iS/™k)Ue[iS/m k ,(i+i)S/m k ) 
satisfies the first one of the two relations above. Furthermore, the second 
one is easily seen to hold for r\ k of this form. □ 

Proof of Theorem {KM and Corollary \2.3X We do this in reverse order, and 
first prove Corollary 12.31 For simplicity of notation we only prove the corol- 
lary for a 2-dimensional Brownian motion, i.e. for the case d = 2. The 
gener al case is the sa me. 

Bv lRootzenl (|l980l ). Theorem 1.2, each marginal process {Hf^ ■ Bj(t), < 



t < T} is tight C([0, T],R), and then also the entire d x d-dimensional 
sequence {H^ ■ Bj(t),0 < t < T,l < ij < d} is tight C([0,T],R dxd ), so 
only stable finite-dimensional convergence remains to be pr oved. We prove 
this in two steps, where the first one follows along the lines of lEootzenl (jl980h 
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and the second step uses the Cramer- Wold device. A final third step uses 
Corollary 12.31 to prove Theorem 2.2. 

Step 1: Let {"0"; i = 1, 2} be adapted processes such that, for i = 1,2, 
(9) sup | / iftds\ -> p 0, 

0<t<T Jo 



and such that 

(10) j\wfds^ p / {^fds, 



t 

2, 







for some ipi, ip 2 > 0, < t < T . To make inverses well defined, we without 
loss of generality can assume that the i/jf(t) are defined also for t > T, and 
such that Equations © and (fTUj) hold with T replaced by S for any S 1 > 0, 
and with ipi{t) = 1 for t > T and i = 1, 2 . This doesn't involve neither 
the result to be proved nor the assumptions, and hence can be done without 
loss of generality. 

Let C[0, oo] = C([0, oo),M) be the space of continuous realvalued func- 
tions defined on [0, oo) and en dowed with t he topology of uniform con- 
vergence on compact sets, see IWhittJ (|1970T ). Let the random variable 



U > satisfy EE/ = 1, and assume the functional / : C[0, oo) — > R is 
bounded and continuous. Further, set r n (f) = f (tpx ) 2 ds + f^ip^^ds, let 
r(t) = lirn n _>. 00 T n (t) = J (^i) 2 ds + f (ip 2 ) 2 ds, and define r" 1 by r^ 1 ^) = 
inf{s : r n (s) > t}. Additionally let TV be a 1-dimensional Brownian motion 
which is independent of T . We first prove that 



/•T"n 1 (-) fTn 1 ^) 

(11) EC//(/ tfdB!+ ^dB 2 )^Ef(W(-)), 

Jo Jo 

for each such U, so that J r " 1( ) iftdBx + J T " 1( ) ip%dB 2 => m W, on C[0,oo). 

Now, define a new probability measure Q by dQ/cflP = C7, and write Eq for 
expectation taken with respect to Q. Then, by Girsanov's theorem (Rogers 
and Williams (2001, Theorem IV 38.5)) there exists an adapted square in- 
tegrable process c = {c\,c 2 ) such that (B(t) = (Bi(t) — J Q c%(s)ds, B 2 {t) — 
Jo c 2 {s)ds) is a Brownian motion under Q. 

Hence, 

(12) EUf( / VidSi + / ^dB 2 ) 

Jo Jo 

^n 1 (-) _ / , T-n 1 (-) 



Eq/( / WdB x + / r/£dB 2 

io Jo 

+ / Vicids+ / ^c 2 ds). 

Jo Jo 



8 CARL LINDBERG, HOLGER ROOTZEN 

Under Q the process same distribution 

as (Rogers and Williams(2001, Theorem IV 34.1)). Further, by Lemma 
12.41 we have that f ipfcids + J^^cids ~~ *p i n C[0, 5 1 ], for any fixed S. 
Since / is bounded and continuous on C[0, oo), these two facts prove (fTT|) . 
and hence mixing convergence on C[0, oo). 

It thus follows from r n -> p r that (r n , J T " 1( ) V™ + /J"" ip%dB 2 )) 
(r, W), and hence, by composing r" 1 with r n , c.f. Billingsley (1999, p. 
145), that 

(13) f VxdB x + f ^dB 2 ^ s W(r(t)) 

Jo Jo 

in C[0,oo), and hence in particular in C[0,T]. 

Step 2: Finite-dimensional stable convergence now follows by standard but 
notationally complicated Cramer- Wold arguments. To lessen complications 
we here only consider two basic cases, and leave the general argument to 
the reader. Thus, first, let ip^(s) = 6jl{ < s <i 4 }-f^{ l j(s) for i = 1,2, with 
< t\ , t 2 < T. Equation ([6]) implies that 

rthti rt/\t2 

r„(t) ^ P r{t) = b\ / Hf^ds + b\ \ H\ 2 ds 
Jo Jo 



so that by ([131) 

rtAtx rtAt2 ptAt\ rtAt 2 

h / Hl x dB x +b 2 \ H{? 2 dB 2 ^ s W{b\ / H^ds+bj / Hl 2 ds). 
Jo Jo Jo Jo 

Now, using elementary properties of Brownian motion together with Rogers 
and Williams(2001, Theorem IV 34.1) we have that W(b{ f* Atl H\ x ds + 

b\ J iAi2 Hl 2 ds) has the same distribution, and the same dependency with 
any J-"-measurable variable, as 

rtAti rtM2 

h / H 1A dW ltl + b 2 / H lj2 dW 1}2 , 
Jo Jo 

for independent Brownian motions W\ } \,Wx t % so that we by (fl~3l) have es- 
tablished that h f* Ml H? 1 dB 1 + b 2 J * A * 2 Hl 2 dB 2 =^ s h f* Atl H hl dW 1A + 
b 2 C t2 H 1>2 dWt, 2 , for any real numbers b\, b 2 . In particular stable 2-dimensional 
convergence of (flft ■ B 1 (t 1 ),H^ 2 ■ B 2 (t 2 )) to (J* 1 H 1A dW 1)U f* 2 H^dW^) 
follows by Cramer- Wold. 

If we instead take V>f = hI{o< s <t 1 }Hf A (s)+b 2 Ii i0 < s < t2 yH^ 1 (s) and ip 2 = 
then, by ©, 

rtAti rtAtiAt2 ftAt 2 

r„(t) ^ P r(t) = bj / Hl 1 ds+2b l b 2 / H 1}1 H 2 , lP { 2 +b 2 2 / H^ds. 

Jo Jo Jo 
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Furthermore, similarly as before, it can be seen that then W(r(-)) has the 
same distribution, and the same dependency with any ^-measurable vari- 
able, as 



/■tAti /-tAti 

( / Hx^dWx,! + / H 1A al 2 dW 2>1 ) 
Jo Jo 

rtM.2 r-tM 2 

+b 2 { H^a^dWx,! + / H 2A al 2 dW 2 ,x) 
Jo Jo 



Reasoning as above we get that 



h [ 1 Hl x dB x + b 2 I 2 H 2X dB\ ^ s 
Jo Jo 

6i I 1 H^a^dW^ + b x [ 1 H 1; xai 2 dW 2jl 
Jo Jo 

+ b 2 #2,102,1^1,1 + b 2 / H 2t x(rl 2 dW- 
Jo Jo 



for independent Brownian motions Wx t x, W%x- Since bx and b 2 are arbitrary, 
this proves stable two-dimensional convergence of (H xl -Bx(tx),H 21 -Bx(t 2 )). 
A general proof of Corollary 12.31 is only notationally more complicated. 
We next use Corollary 12.31 to obtain the conclusion of Theorem 12.21 

Step 3: By Lemma |2~4"1 if satisfies then sup <t<T I Jo ^?j a ids\ — > p 0, 
for all i,j, and hence the general result follows if we can prove that the result 
of the theorem holds for the case when all aj are identically zero. Thus, to 
find the limit of • X{\ one only has to consider 

d 

(^ H ij G 3^ ■ B k }. 
k=l 

However, again by Lemma l2.4( if Hf^ satisfies then 
(14) sup | f H^G hk ds\ ^ p 0. 

0<t<T Jo 

The result then follows from Corollary 12.31 □ 

We now change to a more general setup, from Brownian semimartin- 
gales to general processes (H n ,X n ) which are defined on filtered probability 
spaces * n = (n n , F n ,F n , (Jf) <t<oc)- Here F n is a P n -complete cr-algebra 
and (Ft)o<t<oc is a filtration which satisfies the usual hypotheses (but which 
is not necessarily generated by a Brownian motion) . The following definition 
is key to our goal. We give it for vector valued processes. The definition for 
matrix valued processes is analogous. 

Definition 2.5. Let (X n ) n >x be a sequence of continuous M^-valued semi- 
martingales defined on ^ n , n > 1 and assume that X n =^ X. The sequence 
X n is good if for any sequence of W dxd -valued adapted cadlag stochastic 
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processes (H n ) n >i defined on \I> n such that (H n ,X n ) =>■ (H,X), there exists 
a filtration (Qt) such that X is a semimartingale and H is an adapted cadlag 
process, and {Hfj ■ X™} => {H^j ■ Xj}. 

The following crite r ion is sufficient for goodness, see e. g. Theorem 2.2 



in 



Kurtz and Protterl (jl991al ). 



Definition 2.6. A sequence of continuous M. d -valued semimartingales (X n ) n >i 
is said to have uniformly controlled variations (UCV) if for each n > 1, 
there exist decompositions X n = M n + A n such that 

fT 



supE n {[M n ,M n ] T + I \dA n s \ 

n I JO 



< oo. 



The next theorem combined with Theorem 12.31 will give the asymptotic 
distributions of approximation errors for stochastic integrals. If in addition 
to the conditio ns of the theorem f is boun ded, then the result follows from 
Theorem 3.5 in iKurtz and Protterl (|l991bl ). However, in the present setting 



the result holds also without the boundedness condition, and it is further 
possible to give quite simple proof. In the theorem, = Tq < r ™ < • • • < oo 
are {J^j-stopping times, and r\ n is defined by r] n (t) = t£ < t < t£ +1 . 

Theorem 2.7. Let Y be a continuous M. d -valued {J-t} -semimartingale on 
[0,T], and suppose that f = (/i, . . . , fd) is continuously differentiable. As- 
sume that rj n (t) tends to the identity in probability for t £ [0, T] and let {A n } 
be a positive sequence converging to infinity. Further, set 



U n = \ n J(f(Y)-f(Yo Vn ))dY 

d . 

:= X nY.J (fi(X)-fi(Y°Vn))dYi 



1=1 



and define 

(15) ZfAt) = X n fiYjis) - Yj o ThMdYiis). 

Jo 

Suppose that (Z n ) n >i is good, and that (Z n ,Y) => (Z,Y). Then U n U 
on [0, T], where 



i,5=l 



Since r] n is nondecreasing, pointwise convergence in probability in [0, T], as 
assumed in the theorem, is equivalent to uniform convergence in probability 
in [0, T], Below we will use this without further comment. 



Proof of Theorem 2.1 . For simplicity of exposition, we assume that d = 1. 
By the continuous mapping theorem we have that (Z n ,Y, Y) (Z,Y,Y). 
Since Y is continous and rj n converges uniformly in probability to the unity 
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this in turn can be seen to imply that (Z n , Yorj n , Y) => (Z, Y, Y), e.g. by us- 
ing the Skorokhod translation of convergence in distribution to convergence 
a.s.. 

We now define 

, v f_M - f(y) 

9 ^ V)= —y ' 

where we make the continuous choice g(x, x) = f'(x) when the denominator 
vanishes. The function g is uniformly continuous on [0, T] 2 , so the continuous 
mapping theorem gives that (Z n , g(Y,Y o rj n )) => (Z, f'(Y)). Now, 

U n = \ n J (f(Y) - f(Y o Vn ))dY = J g(Y, Y o r }n )dZ n . 

But since (Z n ) n >\ is good, we have that 

j g(Y, Y o rj n )dZ n =>• j f\Y)dZ, 

which proves the theorem for d = 1. □ 

The next lemma provides a tool for verification of criteria like (jJJ) and © . 
In the lemma we specialize to stopping times (cf. the introduction) defined 
recursively by Tq = and 

(16) r^ +1 = U + A T 



n0(r£), 

for some adapted stochastic process 6. As before let 

(17) Vn(t)=T%, r£<t<rj? +1 , for fc = l,2,..., 

and write £ p = E £ £(s) p ds = £ s p / 2 EB{l) p ds = EB(l) p /(p/2 + 1) so 

that Ei = E Jq B(s)ds = and E 2 = E £ B(s) 2 (is = 1/2. 

In the lemma we will assume that the function a(t); t G [0,T] is locally 
bounded, i.e. that to any e > there exists a localizing stopping time u = u e 
such that a(t Af); t G [0,T] is bounded, and such that F(v < T) < e. In 
particular, if a is continuous on [0,T] then a is locally bounded. 

Lemma 2.8. Assume that a and 9 are adapted processes such that a is 
locally bounded, 8 is strictly positive, and a{t)/6(t) p / 2 is a.s. Riemann inte- 
grable over [0,T], and let and r\ n be defined by dT6j) and (fT7"|) . Set 

oo 

(18) Mt) = n p l 2 Y,a{rZ){B{t) - B(^)) p l {T ^ t<T , +i} , 

k=0 

Further assume that r/ n tends to the identity in probability. Then 

(19) sup | / ip n (s)ds - E p — , /2 ds\ -+ p 0, 
o<t<T Jo Jo 0{s) p/z 

as n —> oo, for p = 1,2. 
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Proof of Lemma \2.8[ If we prove the lemma under the additional restriction 
that a is bounded, then it follows in general, since it then holds for a(t) 
replaced by a(t A v) for any localizing stopping time v, and this in turn 
implies that (I19p holds with probability greater than 1 — e, for arbitrary e. 
Thus we assume in the rest of this proof that a is uniformly bounded, so 
that in particular the expectations exist. 

To ease notation we below sometimes will write Tk instead of tJ} and define 
Tk = T Th . Clearly 

r T k+i fi/ n 9( T k) 
p/ 2 E{ a(T k ){B{t)-B{Tk)Ydt\Fk} = n Pl 2 a{T k ) EB{tfdt 

Jr k JO 



n 



P a ( T k) 



P r,9(T k )P/ 2 + 1 



n> 

Recalling the definition of rj k , 

a(rfc) _ f Tk aori n {s) 
^ v nQ{T k )vl^-^} Q Oo^sy/^ 

and hence 

X k = / lfj n ds - E p / TZT^ds 

Jo JO V°Vn{S) p ' 2 

is a martingale with index set Z+. 

In the following we show that Efc E ((^fc+i - X k ) 2 \F h ) ->• 0- By th e 
functional central limit theorem for martingales (see e.g. Rootzen ( 19831 ). 
Theorem 3.5) this in turn implies that 

f T k f T k a o 77„(s) 

(20) max pCfe = max | / ip n ds - E p / — — -ds\ -+ p 0, 

k k J Q J 0oij n (s)P/ 2 

as n — > oo. Using the Cauchy-Schwarz inequality in the second step, ele- 
mentary properties of Brownian motion in the third, and that (r k +i — T k ) = 
l/{n6(r k )) in the fourth step, we have that 

J2^[(Xk+l-X k ) 2 \f k ] < (r + \ndt) \f k 

k k lA"' 7 * / 

< nPj2 <r k ) 2 {T k+1 - T k ) E[(B(t) - B(T k )) 2p \T k ]dt 

k Jr k 

= — r^n? y~] a(T k ) 2 (T k+ i - Tk) 
p + 1 



\P+2 



E 2 P , a(r k ) ^ a{r k ) 



p + i k " n e( Tk ) P /2+^ k ^ n e( Tk ) P /w 

It follows from the Riemannn integrability of ajQ p l 2 that in the last ex- 
pression above the first factor tends to and that the second tends to 
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a(s)/9(s) p / 2 ds, so that the product tends to zero. This completes the 
proof of pop . 

The assumption that a is bounded and straightforward computation show 
that E J Q T ip^ds is bounded in n, and since furthermore maxfe{r^ +1 — tJ}} — > p 
0, we can apply the Cauchy-Schwarz inequality, to see that 

max sup | / if> n ds\ < I max{r^ +1 — r^} / i^ds I — > p 0, 

k rj}<t<r- +1 JtZ V k JO J 

for n — > oo. Together with (|20p this shows that 

(21) sup | I 4, n ds - E p I Q a °^*l n ds\ -+ p 0. 

0<t<T Jo Jo V°Vn{S) p ' 2 

By assumption a/Q p l 2 is Riemann integrable, and hence 

/ v , /** o,{s) , /"* a o r) n (s) , . 

(22) sup | / - \, 9 ds - / T^ds\ -> a . s . 0. 



o<t<T Jo 0{s) p/2 Jo 0o Vn (s)P/ 2 

The triangle inequality together with (|21 j) and (|22p concludes the proof of 
the lemma. □ 

3. Approximation of stochastic integrals 

We now use the results from the previous section to find the explicit form 
of the asymptotic distribution of the the sum of the errors in approximating 
d stochastic integrals where the integrands are functions of the solution to a 
d-dimensional SDE and where the integrators are the same solutions to the 
SDE. The following condition is used in the theorem. 

Condition 3.1. Let the measurable functions a(-) : R d — > M. d ,f3(-) : R d — > 
R dxd satisfy 

\a{x)\ + \P(x)\ <C(l + |x|), 
where x G M. d for some constant C and 

\a(x) - a(y)\ + \/3(x) - p{y)\ < D\x - y\, 

where x,y G M. d for some constant D. 

This condition ensures that the SDE has an unique continous solution. 
Further, we wi ll need the following lemma, which is given as Lemma 2.5 in 



Rootzenl (11983 ). 



Lemma 3.2. Suppose {Z n } n >i is a sequence of positive discrete time sto- 
chastic processes, adapted to their respective filtrations {J- n } n >i o,nd that r n 
is a stopping time with respect to T n for each n. Then 



J>(Z™| JJU) ^ p 
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implies that 

Theorem 3.3. Let Y be the solution of the SDE 

(23) dY(t) = a(Y(t))dt + (3(Y(t))dB(t), 

where B is a d- dimensional Brownian motion, a,/3 satisfy Condition (|3.ip . 
and Y(0) is independent of B and satisfies EY(0) 2 < oo. Then the error in 
the Euler-type approximation scheme defined by 

U n = n 1 / 2 / f(Y(t)) - f(Y o Vn (t))dY(t), 



where f = (/i, . . . , fd) is continuously differentiable and the grid is given by 
([TBI) with sup tg [ ,T] ^CO < 00 a - s - an d a - s - Niemann integrable, satisfies 

d rt 

/ ^k,m(u)dW k ^ m (u) 
k,m=l ^° 

on [0, T], where 

AkAt) ~m ' 

and W is an d x d- dimensional Brownian motion, independent of B. In 
particular, 

d rt 

sup \U n (t)\ ^ sup | V / A k>m (u)dW k;m (u)\. 

0<t<T 0<t<T k ^ =1 Jo 

Proof of Theorem \3.SX We first assume that the coefficients a and /3 are 
uniformly bounded, and prove that the result holds under this extra as- 
sumption. The general result for unbounded coefficients then follows by an 
easy localization argument which is given at the end of the proof. We again 
write f v instead of F Tv and often suppress the explicit dependence on n and 
e.g. write r v instead of r™. 

Since 1/6 is Riemann integrable, and hence pathwise bounded a.s., and 
su Pie[o,T] Q{t) < oo a.s., i t follow s that rj n tends to t uniformly a.s.. By The- 
orem 5.2.1. in l0ksendall (|2003h there exists a unique t-continuous solution 
Y to Equation ([23D . 

The first part of the proof consists of showing that 

Z n = [\yj(s) ~ Yj o rh( S ))dYi{8) : l<i,j<d 
Jo 

converges jointly with Y. For this we prove that the conditions of Theorem 
I2.2l are satisfied for Hf^ = y/n{Yi—Yiorj n ) and Gj jk = {3j >k - The bounded vari- 
ation part of Yi — Yiorj n can be seen to give contributions which are O p (l/n) 
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and thus, using the triangle inequality and writing l v (s) = 1{ Tv < b <t v+ i}) it 
can be seen that (HI) follows if we show that 



(24) 



n sup 

te[o,T] 



/ l v (s)(3ij(u)dBj(u)ds 



0, 



for 1 < i,j < d. 
Now, 



JO „, Jt„ 



u)ds 



(25) 



Vn yZ l v {s)((3ij(u) - /3ij(T v ))dBi(u)ds 
Jo ~ Jt v 

+ Vn [ V l v {s)pij(T v ){Bi{u) - Bi{T v ))ds. 
Jo „ 



The last term tends to zero in probability by Lemma 12.81 with p = 1, since 
Riemann integr ability of 1/V0 follows from Riemann integrability of 1/6. 

We next show that also the first term on the righthand side is negligible. 
Let C denote a generic deterministic constant whose value may change from 
one appearance to the next. Since t v+ \ is measurable with respect to T v 
it follows from Condition 13.11 Ito's isometry, and the assumption that the 
constants in (1231) are bounded that 



(26) E[ / - ^{r v )) 2 du\T v \ <C E[\Y(u) - Y{r v )\ 2 \T v ]du 

< C / (u- r v )du < C(t v+ i - t v ) 2 . 

J Tv 

Define 

_ rtAr v+ i i-sAt v + 1 

A v (t) = Vn / / (Pi,j(u) - Pij(T v ))dBi(u)ds, 

so that the first term on the righthand side of (|25j) equals A v (t). Using 
Doob's inequality together with the Cauchy-Schwarz inequality in the second 
step and (|26j) in the third step we have that 



E[ sup |A„(t)| \T V ] < ^/n(r v+1 - r v )E[ sup | / - /3ij(T v ))dBi(u)\ \ T v 

T v <t<T v + \ T V <S<T V + 1 J T v 

< C^(t v+1 - Tv )E[ I V+ \Pi,j{u) - p l , j {r v )) 2 du l^} 1 / 2 

J T v 

< C \fn(T v+ i - t v ) 2 . 
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Thus, by the definition (fT6l) . 

J^E[ sup |A„(t)||J-„] < CVriJj 

,, T v <t<T v + l 



T v +1 - T v ) 2 



1 1 

< C\/n—T sup 



n o<t<T V{t) 

According to Lemma [3T2l it follows that ^2 V sup Tv<t<Tv+1 \A v (t)\ —t p 0. Hence, 
sup | VA„(t)| < V sup |A„(t)|-> P 0, 

0<t<T v v T v <t<T v + 1 

which completes the proof that the first term in the righthand side of (|25|) 
tends uniformly to zero in probability. 

Completely similar, but more complex computation show that for any 
indexes i,j,k,l,m, and using Lemma 12.81 with p = 2 for j = m and compu- 
tations similar to (but simpler than) the proof of Lemma 12.81 for j ^ m, 

rt ps ps 

n / ^ / l v (s)l3ij(u)dBj(u) / l v (s)f3^ m (z)dB m (z)f3^ k (s)f3 m ^{s)ds 

JO v Jt v Jt v 

= 71 Pi,j{T v )Pj,k(T v )Pl >m (T v )P m ^(T v ) 
v 

(BjW-BjfaMBnit) - B m (r v )l v (t) + o p (l) 
1 /"* 

^P o / PiAWjAWlAWmAQ/OWMhrn, 
1 JO 

where bj^ m is 1 if j = m and zero otherwise. It follows that also the condition 
© of Theorem O holds, for fl^-(a) = ^{Y^-YiO^s)) and G,-,* = j)k , 
with H itj (s) = fiij/VO and p^ >m = <5j,m, and hence 

(27) {ii-r y j} ME %hk-w h k}, 

l,k=l v w 

since crij,i t m = 1 if j = m, and <7ij,z, m = otherwise. 

Arguments similar to those above show that {-f^,- ■ 5^} has uniformly 
controlled variations and hence are good. Stable convergence implies that 
the left hand side of Equation [57] converges jointly with Y. The first con- 
clusion of the theorem now follows from Theorem 12.71 f° r the case when the 
coefficients are bounded. 

To remove the restriction that the coefficients are bounded, for general 
ai,/3i,j define c oefficients a\ = on A c and (3f j = faj A c. Theorem 5.2.1. in 
0ksendal (120031 1 still yields unique f-continuous solution Y c to Equation (|23p 



for these functions. Let JJ n,c be defined from a?, • in the same way as U n 
is defined from «j, With obvious notation, we have already proved that 
U n ' c => s U c , as n -> oo for each fixed c. Since P(sup t6 [ 0)T ] |l" c (t) - Y(t)\ > 
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0) — > 0, as c — > oo also U c =^ s U. Further, 
limsupP( sup \U n ' c -U n \> 0) 

n te[0,T] 

< P(inf{f : max{max{a i (lt)} ) max{A )J -(lt)}} > c} < T) — > 0, 
as c — ?• oo. Hence, Theorem 3.2 in iBillingslevI (|l999l ) gives that C/ n [7, 



which proves that the first result of the theorem holds also for the general 
case. 

The second conclusion follows from from the first by the continuous map- 
ping theorem, since the supremum mapping is continuous. □ 

4. Designing the error in approximations of stochastic 

integrals 

In deciding on which approximation scheme to use to compute a stochastic 
integral - or, to decide on a hedging strategy - one has to balance the error 
with the number of intervention times N = N n = maxjfc; t£ < T}. In this 
section we will investigate two such schemes. The first one could be called 
the "no bad days" strategy, and simply consists in choosing the stopping 
times {t/%} where the stochastic integral is evaluated - or the times when 
the portfolio is rehedged - in such a way that the error is a Wiener process. 
In the second strategy we bound the expected number of evaluation times 
and minimize the asymptotic standard deviation of the approximation error 
under this restriction. 

The setting of this section is the following: Suppressing the superscript n 
the stopping times are given by (fTBj) . i.e. tq = and 



(28) ™ = ( Tk+ ^k)) AT 

with 9 adapted and positive, and the distribution of the approximation error 
e(t) satisfies 

(29) V^(*)=> fJ^LdWit), 

Jo ^6{s) 

for some adapted process f(s)>0 and Wiener process W which is indepen- 
dent of and /. Here it should be noted that (|29p is more general than it 
looks at first. E.g. the approximation error in Theorem 13.31 satisfies this for 

It is straightforward to find the asymptotic number of evaluation times. 

Proposition 4.1. Suppose that 9 is Riemann integrable a.s. and that 
info<t<T 9{t) > a.s.. Then 

rT 



lim — - = / 9(t)dt, a.s. 

n->oo n Jq 
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//, in addition, E[sup 0<t<T 9(t)] < oo, then 



n, r T 

lim E— = / E9(t)dt. 



Proof of Proposition \4-l\ Suppose first 9 is of the form 

k 

(30) 0(t) = ^l[ Oi ,a i+1 )(*), 



i=0 



for some random variables 9{ > and = ao < a\ < ... < = T, and 
with l[ 0ij0i+1 ) the indicator function of the interval [aj,aj+i). For each u, it 
is easily seen that the number of intervention times in the interval [a^, aj+i) 
is n9i(ai + i — ai) + 0(1), and hence 



N 1 1 /" 

-2 = Y j 9 i {a i+l -a l ) + 0(-)= / e(t)£ft + 0(-)-». / 



T 

9{t)dt, 



as n — t- oo. If # < and is of the form (|30p then, with obvious notation, 
N n (9) < N n (9)+0(l), and the corresponding bound with all the inequalities 
reversed is also true. 

Now, by assumption 9 is Riemann integrable, and hence can be approxi- 
mated arbitrarily well from below and above by functions of the form (|30p . 
This proves the first assertion of the proposition. 

Furthermore, N n /n < T sup 0<t<T 9{t) + 1/n, and hence the second asser- 
tion follows from the first one by dominated convergence. □ 

In the rest of this section we assume that we "are in the asymptotic 
regime", i.e. that n is so large that we, to the degree of approximation 
needed, may assume that the limits above can be replaced by equalities. 
Thus, below we will assume that 



(31) EN = nl E0(/)r//. eft) = 4= I ^= dW(t), 

In 



T 1 rt 

o 

so that in particular Ee(t) 2 = ^ J * ds. 

The no bad days strategy: It is at once seen, supposing that f 2 is Riemann 
integrable, that if we choose 9{t) = cf(t) 2 , for some constant c, then 

e{t) = -}=W(t) 
yen 

and 



r T 

EN = cn / Ef 2 (s)ds. 
Jo 
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Thus, in a financial setting, with this choice of there are no "days" where 
the hedging error grows quicker than during other days, and hence a trader 
can sleep equally well (or equally badly!) each night. 

Minimal standard deviation: We will now, supposing that / is Riemann 
integrable, show that the solution of the optimization problem 



inf {yjEe 2 {T) : EN < C} 

{0-. e>o, adapted} 

is given by 6(t) = Cf(t)/(n Jo Ef(s)ds). For this choice 



/ fcEfds , 
EN = C, e(t) = y J ° c / V /'/"'• 



o 



Thus in particular, for the optimal strategy the standard deviation is y / Ee(T) 2 
f^Efds/VC. 

Now, write 9 = n9. With this notation Ee(T) 2 = E J Q T f 2 /6ds and the 

restriction is E J Q T 6 < C. Applying the Cauchy-Schwarz inequality twice, it 
follows that 

(e£ fdsf < (e^J £ p/ed S] j £ eds 

< E( I f/eds)E( [ 6ds), 
Jo Jo 

and hence 

, N2 {E^fdsf 

However, above we have seen that § = ne = C f /($^ Efds) achieves this 
bound, and hence is the optimal choice. 

5. Application to hedging 

An important application of the results in the previous section is to hedg- 
ing of financial derivatives. Here we treat the simplest Black-Scholes model 
and only give a brief comment on more complicated problems. The limit 
distribution of the Black-Sch oles hedging error for e quid istant deterministic 
grids h as been studied in e.g. iBertsimas et al. (l200(]h and lHayashi and Myklandl 



(2005 ). (We have no t been able to follow the proof of Theorem l.b in 



Bertsimas et al.l (120001) - specifically we c ould not understand the use of 



Lemma 5.1 from lDuffie and Protterl (119921 ) 



We distinguish between complete and incomplete financial markets. In 
complete markets, all derivatives can be replicated (hedged) perfectly by 
trading in a self-financing way in the underlying and a money market ac- 
count. The approximation error distribution we analyze is here the total 
hedging error. In an incomplete market, an investor who hedges a contract 
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will still choose a hedging portfolio which is, in some sense, optimal for her 
purposes. In this case, the error we obtain is relative to this optimal hedging 
portfolio. We give now an application of the results in the previous section 
to hedging in the complete Black-Scholes market. 

We assume that a stock S follows the Black-Scholes model. In other 
words, we model the stock as a geometric Brownian motion, which has the 
dynamics 

dS(t) = nS(t)dt + aS(t)dB(t), 

for //, a > 0, where B is a Brownian motion, and S(0) = s > 0. Further, we 
have a risk-free money market account with dynamics 

dR(t) = rR(t)dt, 

for r > 0, where -R(O) = 1. It is well known that the price of a so-called call 
option with payoff max(S(T) — K, 0) at the deterministic terminal time T, 
for some strike price K, is at time t 

n(t) = *(d+)5(t) - Ke'^^id-) 

where <3? denotes the standard normal cumulative distribution function and 



d±(t) 

Now, if we set 



Iog^ + (r±4)(T-t) 



Y(t) 



Sit) 
R(t) 



and / = (*(</+), -$>(d-)Ke- rT ), we get that 

n(t) = J f(Y(t))dY(t) 

gives the self-financing price process of the call option. This is of the form 
considered in Theorem 13.31 with d = 2 and = crS(t), and all other 

/3-s equal to zero. Thus, using the stopping times (fTB]) . Theorem 13.31 gives 
that the hedging error satisfies 

Vn(n(t)-n OJh (f)) ^ s y ^-(s)a 2 S(s) 2 /^20(s)dW( S ) 

= r pg^w, 

7o V 2 ^( S )( T - S ) 

with = d<&(t)/dt the standard normal density function. 

Consider now an investor who hedges a call option, but who only adjusts 
her hedge at some stopping times {rk}k>i of her own choosing. If she want to 
have a " uniform" increase of the error and make it approximately a Brownian 
motion, she should use the "no bad days" strategy from the previous section. 
This would mean that she would use the stopping times (|28h with 9{t) = 
c4>(d + (t)) 2 a 2 S(t) 2 /(2(T — t)). However, this leads to a (purely) technical 
difficulty: B(t) tends to as t -> T if S(T) G E \ K and to oo if 5(T) = A". 
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This means that the assumption of a.s. Riemann integrability of 1/9 is not 
satisfied on [0, T], nor is the assumption that sup t6 r T ] 8{t) < oo on [0, T]. A 
theoretical (and in fact also practical) solution is to instead only evaluate the 
hedging strategy up to a constant time V < T, with V close to T. Theorem 
(|3.3p gives that the hedging error up until V for large cn then approximately 
is distributed as W{t)/y/cn. 

Alternatively, the minimum standard deviation strategy and the same 
reasoning as above leads to choosing 

v^fr-f) h v '2(r - 5) 1 

where C is the expected number of evaluation times. This yields the ap- 
proximate distribution 

for the hedging error, for n large. 

It is now completely straightforward to add one or more stocks to the 
portfolio and, using e.g. that Jj fi/VOdWi + jj f 2 /y/6dW 2 has the same 
distribution as J* * \//i + /i / V@dW, to find the optimal stopping times and 
the resulting error when the hedges for all of the stocks are adjusted at the 
same time points. This is how portfolio hedging is done in practice. We 
leave these calculations to the reader. 



References 

Aldous, D. J., Eagleson, G. K. (1978): On mixing and stability of limit theorems. 

Annals of Probability 6, 325-331. 
Bertsimas, D., Kogan, L., Lo, A. W. (2000): When is time continuous? Journal of 

Financial Economics 55, 173-204. 
Billingsley, P. (1999): Convergence of Probability Measures. John Wiley & Sons, Inc., 

New York. 

Broden, M., WlKTORSSON, M. (2010): Hedging errors induced by discrete trading under 
an adaptive trading strategy. Preprint. 

Duffie, D., Protter, P. (1992): From discrete to continuous time finance: weak con- 
vergence of the financial gain process. Mathematical Finance 2, 1-15. 

Fukasawa, M. (2011) Discretization error of stochastic integrals. To appear in Annals of 
Applied Probability. 

Geiss, S., Toivola, A. (2009): Weak convergence of error processes in discretizations of 

stochastic integrals and Besov spaces. Bernoulli 15, 925-954. 
Gobet, E., Temam, E. (2001): Discrete time hedging errors for options with irregular 

payoffs. Finance and Stochastics 5, 357-367. 



22 



CARL LINDBERG, HOLGER ROOTZEN 



Hayashi, T., Mykland, P. A. (2005): Evaluating hedging errors: an asymptotic ap- 
proach. Mathematical Finance 15, 309-343. 

Jacod, J., Protter, P. (1998): Asymptotic error distributions for the Euler method for 
stochastic differential equations. Annals of Probability 26, 267-307. 

Kallenberg, O. (2002): Foundations of modern probability. Springer Verlag, New York. 

Kurtz, T. G., Protter, P. (1991): Weak limit theorems for stochastic integrals and 
stochastic differential equations. Annals of Probability 19, 1035-1070. 

Kurtz, T. G., Protter, P. (1991): Wong-Zakai corrections, random evolutions, and 
numerical schemes for SDEs. In Stochastic Analysis 331-346. Academic Press, New 
York. 

Kurtz, T. G., Protter, P. (1996): Weak convergence of stochastic integrals and differ- 
ential equations. 1995 CIME School in Probability. Lecture Notes in Math. 1627, 1-41. 
Springer Verlag, Berlin. 

ROGERS, L.C.G. and Williams, D. (2001) Diffusions, Markov Processes and martin- 
gales, Vol 2, Ito Calculus. University Press, Cambridge. 

Rootzen, H. (1983): Central limit theorems for martingales via random change of time. 
In Essays in honour of Carl Gustav Essen, Eds. Hoist and Gut, Uppsala University, 154 
- 190 (available as item [21] on http://www.math.chalmers.se/~rootzen/lista.html). 

Rootzen, H. (1980): Limit distributions for the error in approximations of stochastic 
integrals. Annals of Probability 8, 241-251. 

Tankov, P., Voltchkova, E. (2009): Asymptotic analysis of hedging errors in models 
with jumps. Stochastic Processes and their Applications 119, 2004-2027. 

Whitt, W. (1970): Weak convergence of probability measures on the space C[0, oo). 
Annals of Mathematical Statistics 41, 939-944. 

0KSENDAL, B.(2003): Stochastic differential equations: an introduction with applications, 
6th ed. Springer Verlag, Berlin Heidelberg. 

C. Lindberg: Mathematical Sciences, Chalmers University of Technology 
and Gothenburg University, SE-412 96 Goteborg, Sweden 

H. Rootzen: Mathematical Sciences, Chalmers University of Technology 
and Gothenburg University, SE-412 96 Goteborg, Sweden; 

|HTTP : / / WWW. MATH . CHALMERS . Se/~ROOTZEN/ | 



